An Isocurvature Mechanism for Structure Formation 
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We examine a novel mechanism for structure formation involving initial number density fluctuations 
between relativistic species, one of which then undergoes a temporary downward variation in its 
equation of state and generates superhorizon-scale density fluctuations. Isocurvature decaying dark 
matter models (iDDM) provide concrete examples. This mechanism solves the phenomenological 
problems of traditional isocurvature models, allowing iDDM models to fit the current CMB and 
large-scale structure data, while still providing novel behavior. We characterize the decaying dark 
matter and its decay products as a single component of "generalized dark matter". This simplifies 
calculations in decaying dark matter models and others that utilize this mechanism for structure 
formation. 



Conventional models for structure formation utilize ei- 
ther initial density fluctuations that then grow by grav- 
itational instability or initial stress fluctuations between 
the matter and radiation which push the matter into 
gravitationally unstable configurations. We call these 
adiabatic and traditional isocurvature models, respec- 
tively. The possibility remains that the origin of struc- 
ture lies elsewhere. We consider here an origin from a 
temporary change in the equation of state of the dark 
matter and give a concrete example in the form of an 
isocurvature decaying dark matter (iDDM) model for 
structure formation. 

Traditional Isocurvature Models. — In traditional 
isocurvature models, the initial density balance comes 
at the expense of number density or "entropy" fluctu- 
ations between the matter and radiation. These mod- 
els suffer from various problems that stem from the 
lack of a natural timescale besides the expansion time 
in the radiation-dominated early universe. Phcnomcno- 
logically, models with scale-invariant entropy fluctua- 
tions differ from their adiabatic counterparts with scale- 
invariant curvature fluctuations through: (1) an under- 
production of large-scale structure relative to large-angle 
CMB anisotropies jij and (2) diminished acoustic peaks 
appearing at small angles Q. The former is in conflict 
with the observed large-scale structure, once normalized 
to large-scale anisotropies, and the latter is in conflict 
with recent measurements of degree-scale anisotropies. 

Variants of the basic isocurvature model have been 
proposed to alleviate these problems. Peebles Q intro- 
duced a blue-tilt to the entropy power spectrum to ad- 
dress the lack of large-scale structure and the falling spec- 
trum of CMB anisotropies. Unfortunately, such models 
are in moderate conflict with the slope of the COBE 
anisotropy spectrum Q|. The basic problem with the 
spectrum is one of timescale. This model forms struc- 
ture through the residual stress perturbation remaining 
when the density fluctuations in a fluid with relativistic 
pressure p r = p r /3 are balanced off those in a pressure- 
less fluid. These stresses move matter around until the 
perturbation crosses the horizon or the radiation den- 
sity becomes negligible. For wavelengths smaller than 



the horizon at matter-radiation equality, the amount of 
time the stresses have to act is a decreasing function of 
wavenumber k, leading to a falling spectrum of CMB 
anisotropies and large-scale structure. 

Cosmological defects provide an alternative means of 
generating structure from isocurvature initial conditions 
that reduce these problems of spectral shape. These 
models balance seed fluctuations in the defects off or- 
dinary matter to establish isocurvature initial condi- 
tions. Unlike the matter-radiation isocurvature models, 
the temporal behavior of the stresses scale with the hori- 
zon crossing time || . Thus the effect of the stresses on 
the ordinary matter is the same for all scales. The other 
problems of isocurvature models remain. For example, 
these models tend to underpredict large-scale structure 
(see for recent assessments). Part of this problem 
would seem to be common to all isocurvature models. 
It comes from the relationship between the Sachs- Wolfe 
temperature anisotropy AT/T and the Newtonian curva- 
ture <&, AT/T « — 2$, which should be compared with 
AT/T w <E>/3 in adiabatic models. Defect models fare 
even worse since their vector and tensor modes provide 
additional sources of anisotropy. Furthermore, the super- 
horizon growth of the curvature leads to forced acoustic 
oscillations whose features are shifted to smaller angles 
relative to adiabatic models. This would also seem to be 
a generic feature of isocurvature models since the curva- 
ture can only grow from its vanishing initial condition. 
The observed rise of the anisotropy power spectrum at 
degree scales would be difficult to explain in defect or 
other similar isocurvature variants 

Novel Isocurvature Mechanism. — Nonetheless, these 
problems are not fundamental to isocurvature initial con- 
ditions but rather to the choice of the stress perturbation 
history that generates the structure. In both examples 
above, the problem is that horizon crossing for the per- 
turbations sets the timescale over which the stresses act. 
If the stresses on all superhorizon scales could be turned 
on and then off uniformly, then the universe would be 
left with constant scale-invariant curvature fluctuations 
similar to an adiabatic model. 

A dark matter species that undergoes a variation in its 
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equation of state provides such a mechanism. Consider 
the case where a dark matter species goes non-relativistic 
and then decays: the equation of state for this matter and 
its decay products begins at u> d dm = Pddm/Pddm = 1/3, 
dips toward zero, and returns to 1/3 (see Fig. |l| lower 
panel). Another important aspect of this model is that 
since the perturbations are in a species that is originally 
ultra-relativistic, isocurvature conditions require balanc- 
ing perturbations from the other relativistic species, in- 
cluding the photons in the CMB. Because the balance 
is through species with the same equation of state ini- 
tially, vanishing density perturbations imply vanishing 
stress perturbations as well. The initial stress-energy 
tensor of the total matter is completely homogeneous and 
isotropic. 

As the dark matter becomes non-relativistic, the ini- 
tially counterbalancing stress fluctuations become unbal- 
anced as Sp ddm drops below Sp ddm /3. They move mat- 
ter around and form density or curvature fluctuations. 
When the dark matter then decays back into radiation, 
the stresses regain their balance and stop forming curva- 
ture fluctuations. The process that generates curvature 
in this model in fact has an exact solution [^j in the sim- 
ple case that the fluctuations in the decaying dark matter 
are balanced off a single radiation component (r). The 
curvature in the comoving gauge £ is directly related to 
the stresses JlIJ and is given by 

tl , cr p ddm + p ddm a 

Cw = o ; +C(0), 



3 p + p 
where a is the scale factor and 
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Isocurvature models have £(0) = by definition. The 
generation of curvature is most effective if the DDM 
species comes to dominate the total energy density be- 
fore the decay. The curvature remains constant after the 
decay in all cases since the decay products and ordinary 
radiation redshift in the same way leaving p ddm / p con- 
stant until other matter species become important. 

This mechanism solves the problems isocurvature mod- 
els have with the acoustic peaks in the CMB. The con- 
stancy of the potential outside the horizon during radi- 
ation domination insures acoustic phenomenology that 
is similar to the adiabatic model: scale-invariant oscilla- 
tions before diffusion damping and a cosine series of har- 
monic peaks. As noted by many authors in the adiabatic 
variant of the DDM scenario (e.g. P,|i"o||), the decay into 
relativistic species also delays matter-radiation equality 
and allows a high-O m model to look like a low-O m model 
with respect to the large-scale structure power spectrum 
and CMB anisotropies. 

Unfortunately, this mechanism does not automatically 
solve the other problem of isocurvature models: the ra- 
tio of large-scale anisotropies to large-scale structure. In 
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FIG. 1. Time evolution of the comoving curvature £, 
Newtonian curvature $, effective temperature perturbation 
AT/T, and equation of state uiddm in m = 5 keV, r = 5yrs, 
flm — 1, h — 0.5, Qth 2 = 0.02 universe. 



fact, it can actually make the agreement worse. If the 
DDM fluctuations are balanced by CMB fluctuations ini- 
tially, then the CMB temperature will be lower in re- 
gions where the DDM density is higher. Since these 
are the sites of potential wells when the DDM become 
non-relativistic, those cold photons will lose more energy 
climbing out of the potential wells after last scattering 
and become colder. The complete gravitational redshift 
effect is 
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The solution is clear. The initial conditions must involve 
the photons and DDM together compensating the other 
species of radiation. The photons will then be initially 
hotter at the sites that form potential wells. 

iDDM Models. — For definiteness, let us now examine 
the phenomenology of the iDDM class of models. The 
important aspect is that the DDM scenario allows the 
required variation in the equation of state. To highlight 
this property, we model the decaying particle and its de- 
cay products as a single component of "generalized" dark 
matter of the type introduced in |l2|; we label it DDM 
here. The DDM is described by an equation of state 
Wddm, a sound speed Cddm, and a viscosity parameter 
addm- The equation of state of the DDM is given implic- 
itly in terms of the mass m and lifetime r of the particle, 
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w 1 = -[l + (a/a ni )^}- 1 ^, (5) 

with p = 0.872 and a nr = 8.3 x lO^^m/keV)" 1 ^/ 7 ^)- 
Here T/T 7 is the ratio of temperatures of the DDM and 
photons while the DDM was relativistic. We assume here 
that the DDM accounts for one of the usual neutrino 
species so that T/T 7 = (4/11) 1 / 3 . The sound speed is 
given by c^ dm = w ddm - (dw ddm /dlna)/(3 + 3w ddm ) and 
the viscosity parameter a ddm = Wddm following (l^] . We 
have checked this approximation against explicit calcula- 
tions of an adiabatic DDM model with the techniques of 
^ . The approximation has the practical benefit of being 
simple to implement (c.f. |]l3|) and the pedagogical value 
of highlighting the important aspects of this mechanism. 

The critical parameter is m 2 r pJ], since it determines 
the fractional increase in the total radiation density due 
to the decay p r — ► p r (l + 0.15r7i^yTy/ 3 ); as we have seen 
changes in the energy density of the DDM relative to the 
other components is what drives curvature generation. 

We establish isocurvature initial conditions of the form 
Spddm _ $Pj_ _ 3 5p m 

Pddm P-y 4 p m ' 

8p„ = -(Spddm + <Vddm + Sp m ) , (6) 

where p v is the density in the remaining two neutrino 
species and p m is the density in baryons and CDM with 
all other metric and matter perturbations zero. More 
complicated initial conditions can alter the ratio between 
AT/T and Despite these simple initial conditions, 
some care must be taken to insure numerical stability in 
the evolution. In the synchronous gauge, we recommend 
the use of the variables Ht and — Hl — Ht/3 in the 
notation of pj ]. 

The time evolution of superhorizon scale perturbations 
in an example model with m = 5keV and r = 5 yrs is 
shown in Fig. [l] (upper panel). As discussed above, the 
comoving curvature £ grows rapidly as the equation of 
state Wddm dips below its relativistic value. On the other 
hand, the curvature in the Newtonian gauge <I> is finite 
and constant before the decay. In traditional isocurva- 
ture and adiabatic models, the £ and $ are simply pro- 
portional. The difference here is that the radiation pos- 
sesses substantial density perturbations which leads to 
anisotropic stress perturbations in the neutrinos of order 
7r = {kr/) 2 0(6p ddm / p ddm ) where 77 is the conformal time; 
balancing anisotropies in the photons are prevented by 
Compton scattering. The relativistic Poisson equation 
involves curvature sources of order ir/(ki]) 2 leading to 
constant Newtonian curvature perturbations. However 
once the dark matter decays, the Newtonian and comov- 
ing curvatures again become proportional to each other. 

The second interesting feature is that, unlike adiabatic 
models, the comoving curvature £ grows at the second 
transition to matter domination a cq . The initial condi- 
tions (||) implies an entropy fluctuation between the non- 
relativistic matter p m and the combined radiation after 



the decay. This entropy fluctuation causes stress fluctu- 
ations when the non-relativistic matter comes to domi- 
nate the universe. This behavior is identical to matter- 
radiation isocurvature models and implies that the phe- 
nomenology of the complete model will be a combination 
of traditional isocurvature and adiabatic models. Finally, 
as desired, the initial temperature perturbation AT/T 
cancels part of the redshift effect from the growing New- 
tonian potential leaving smaller CMB anisotropies for a 
given matter density fluctuation. 

The CMB anisotropies of models with m = 5keV and 
r =3, 5, 8 yrs are shown in Fig. |^ and compared with 
the current data. Models with the same to 2 t are approxi- 
mately degenerate. The large angle anisotropies have the 
spectral shape of traditional isocurvature models with a 
decline toward smaller angles due to the matter-radiation 
transition. However, the acoustic peaks show a dis- 
tinctly adiabatic pattern with peak positions in a classic 
1:2:3... series and high odd peaks. The result is a dip 
in the anisotropies at intermediate scales that is an inter- 
esting signature of such models. The CMB polarization 
however mimics adiabatic models predictions. 

The power spectrum of the large scale structure in the 
same models are given in Fig. [| and compared with the 
data |Q. The shape of the power spectrum matches 
the data adequately; for example, the ratio of power at 
the 50/i _1 Mpc to 8/i _1 Mpc scale a 50 /a s = 0.16 in 
the t = 3yrs model. Relative to an adiabatic model 
with the same parameters, these have more large com- 
pared with small scale power due to growth at the second 
matter-radiation transition. Note that the observations 
may be slid up and down to account for an unknown 
galaxy bias. On the other hand, the model curves are 
COBE-normalized through the fitting form of jl6| and 
predict erg = 0.62 for the 3yrs model. This value would 
reproduce the present-day cluster abundance adequately 
(e.g. and at high-redshift marginally (e.g. (lljl ). The 
high value of Q m in these examples is observationally dis- 
favored by recent determinations of the the luminosity 
distance to high-redshift supernovae (|^j2l[] ; if these pre- 
liminary indications are borne out by future measures, 
lower Q m variants of these models can be considered. 

Discussion. — Perhaps more interesting than the de- 
tails of any specific model of this type is the lesson 
their mere existance teaches us. That the current CMB 
and large-scale structure data is consistent with a high- 
density model with scale-invariant initial isocurvature 
perturbations warns us against overinterpreting the cur- 
rent data. It also suggests that if future observations 
reveal phenomenology that is close to but not precisely 
predicted by standard adiabatic models, we should not 
necessarily abandon the search for new paradigms for 
structure formation. Our example shows that the acous- 
tic peak locations do not actually discriminate between 
adiabatic and isocurvature initial conditions in general. 
They do however relate to the role of causality in the 
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FIG. 2. CMB anisotropies for models m = 5keV and 
r =3,5,8 yrs (short-dashed, solid, long-dashed lines) in an 
On = 1, h = 0.5 and Q b h 2 = 0.02 universe compared with 
the current CMB data jnj with la error boxes shaded ac- 
cording to area. 
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FIG. 3. Matter power spectrum for the same models as 
Fig. |^ compared with galaxy survey data from [Q with la 
error boxes. Note that the data may be shifted up or down 
to account for the unknown galaxy bias. 



generation of perturbations as discussed in To be 

explained causally, the scale-invariant isocurvature "ini- 
tial" conditions employed here still require a period of 
superluminal expansion, like that provided by inflation. 
The models considered in the last section require that 
fluctuations in the chemical potential of say the e and fi 
neutrinos balance temperature fluctuations and may be 
difficult to arrange in inflationary models. However the 
general mechanism is of wider interest. Exploration of 
alternatives such as these helps to sharpen the questions 
that can be asked of the data. 

In summary, the proposed isocurvature mechanism for 
structure formation solves the problems with the spectra 
and relative amplitude of the matter and CMB fluctu- 
ations associated with traditional isocurvature models. 
The key aspects of this mechanism are that (1) joint 
fluctuations in the energy density of some initially rel- 
ativistic component and the photons are balanced by the 
remaining radiation and (2) the equation of state of the 
DDM drops and then returns to the ultrarelativistic value 
of Wddm = 1/3 before the scales relevant to large-scale 
structure and CMB anisotropies enter the horizon. We 
have shown that this situation can be realized with a dark 
matter particle in the keV mass range that decays on the 
timescale of a year. Such models are in agreement with 
the current data but have novel features that are testable 
with the upcoming generation of experiments. Their exis- 
tance calls into question widely-held beliefs about isocur- 
vature models and cautions agains overinterpretation of 
current and future data sets. 
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